Wigner-PDC description of photon entanglement as a local-realistic theory 
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Regardless of whether the zeropoint field (ZPF) of electromagnetic vacuum fluctuations is con- 
sidered as real or as merely virtual, here we show that the Wigner description of PDC-generated 
(Parametric Down Conversion) photon-entanglement (developed in 0-0]) can be formulated as en- 
tirely local-realistic. This involves reinterpreting the expressions for the detection probabilities, by 
means of an additional mathematical manipulation; though such manipulation seemingly provides 
enough freedom to guarantee consistency with the expectable, experimentally testable behavior of 
detectors, this is, in any case, irrelevant in relation to our main result, of a purely mathematical 
nature. This said, of course the PDC model is by construction restricted to a certain subset of QED- 
states, obtained directly from a mix of the vacuum state (hence one with positive Wigner function) 
with a quasi-classical signal (the laser), and a time-evolution governed by a quadratic Hamiltonian 
(hence one that preserves the positivity of the Wigner function): it is only for that subset that the 
local-realistic analogy can be implemented. Some additional questions are also addressed; amongst 
them, we provide an extremely simple interpretation on that apparently awkward "subtraction" of 
the average ZPF intensity at the Wigner image of the detection process, an element stemming from 
Glauber's purely quantum mechanical expression. On the other hand, the immediate connection 
of this model with the phenomenon of "detection probability enhancement", with the potential to 
reconcile local realism with the experimental violations of inhomogeneous Bell inequalities, cannot 
not be considered anything short of strong evidence of its physical relevance. 

PACS numbers: 03.65.Ud, 03.67.Mn, 42.25, 42.50.Xa 



I. INTRODUCTION 



The truth is that so far no conclusive proof of the vio- 
lation of local-realism has ever been obtained (in this re- 
gard, the detection loophole is also present in the latest 
important development [l|); perhaps it is time to con- 
sider the possibility that not all states allowed by QM 
may have a physical counterpart, in particular those that 
would yield correlations defying local-realism. Perhaps 
is also time to start giving credit to models that, even 
though being quantum-mechanical, may include that re- 
striction built within their formulation; here we propose 
a step in this direction. This step is necessary, but of 
course not conclusive: in coherence with what we just 
said, what we treat here is just that, a "model". We 
hope, nevertheless, that it serves the purpose of encour- 
aging further work on the issue. 

The Wigner picture of Quantum Optics of photon- 
entanglement generated by Parametric Down Conversion 
[1| was developed some years ago in a series of papers [1- 
|9[; recently, the approach has been revitalized produc- 
ing another stream of very interesting results (lol . ITlj . 
It departs from an stochastic electrodynamical descrip- 
tion (hence, based on continuous variables: electromag- 
netic fields defined at each and every point of space), 
and by use of the so-called Wigner transformation [12] it 
acquires a form where all expectation values depend on 
a certain probability distribution for the value of those 
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fields in the vacuum. Such a picture clearly differs, on the 
other hand, from the usual one in Quantum Information 
(QInf), based on a discrete description of the photon (a 
particle) and not on a set of fields in a continuum space. 
For instance, in a polarizing beam splitter (PBS), empty 
polarization channels at either of the exits of the device 
are filled with the random components of the vacuum 
Zero-Point field (ZPF); these random components can 
for instance give rise to the enhancement of the detec- 
tion probability for certain realizations of the state of the 
fields (certain photons, in the QED language): see note 
. Phenomena such as this last come as a big difference 
with the customary model of the set polarizer-detector is 
treated just as a "black box" , able to extract polarization 
information in principle without the (explicit) interven- 
tion of any additional noise. 

The Wigner-PDC framework provides an alternative, 
apparently local-realistic explanation of the results of 
many typically quantum experiments, the result of a 
measurement depending on (and only on) the set of hid- 
den variables (HV) inside its light cone: in this case, 
both the signal propagating from the source and the ad- 
ditional noise introduced by the ZPF at intermediate de- 
vices and detectors. Within the first series of papers 
those experiments included: frustrated two photon 
creation via interference Q ; induced coherence and indis- 
tinguishableness in two-photon interference [1, Q ; Rarity 
and Tapster's 1990 experiment with phase-momentum 
entanglement Franson's (original, 1989) experiment 
quantum dispersion cancellation and Kwiat, Stein- 
berg and Chiao's quantum eraser From the most 
recent one [13, [HI we can also add, on one side, amongst 
quantum cryptography experiments based on PDC: two- 
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qubit entanglement and cryptography [l^ and quantum 
key distribution and eavesdropping [l^l; on the other, an 
interpretation for the experimental partial measurement 
of the Bell states generated from a single degree of free- 
dom (polarization) in [ll]. We have also had access to 
some promising (but still unpublished) work on hyper- 
entanglement [14]. 

Perhaps needless to say, this "local realistic" picture 
of the quantum experiments takes place within the lim- 
its (detection efficiencies) where it is already well ac- 
knowledged that such an explanation may exist; more- 
over, those limits arise as natural consequences of the 
theory: they are simply limiting values of detection rates 
for states of light that are ultimately generated from the 
vacuum. This last feature is what makes the Wigner ap- 
proach so interesting; on the other hand, above we have 
used "apparently" because that local-realistic interpre- 
tation, even within the corresponding limits of detection 
rates, was until now not devoid of difficulties, in particu- 
lar related to the detection model (so far merely a one-to- 
one counterpart with Glauber's original expressions [l5|): 
as a result of the normal order of operators, there average 
intensity due to vacuum fluctuations is "subtracted" , a 
subtraction that seems to introduce problems related to 
the appearance of what could be interpreted as "negative 
probabilities". This is the issue that interests us here, 
one that, as expectable, motivated the proposal of sev- 
eral modifications upon the expressions for the detection 
probabilities: for instance, as early as in Q, "our the- 
ory is also in almost perfect one-to-one correspondence 
with the standard Hilbert-space theory, the only differ- 
ence being the modification in the detection probabil- 
ity that we proposed in relation...". Those modifications 
ranged from the mere inclusion of temporal and spatial 
integration Q to the proposal of much more complicated 
functional dependencies all of them seeming to 

pose their own problems; for instance in [l6j . a departure 
from the quantum predictions at low or high intensities. 

Our route here is a different one however: we do not 
propose any modification of the initial expressions for 
the detection probabilities (in one-to-one correspondence 
with the initial quantum electrodynamical model), but 
just explore the possibility of performing some conve- 
nient mathematical manipulation that casts them in a 
form consistent with the axioms of probability, hence one 
consistent with local-realism. Following for instance @ 
(see also [H), quantum mechanical detection probabili- 
ties, single and joint, can respectively be expressed as 

Pi cx {li - Io,i) 

= I {h{a,a*)-Io.{)W{a,a*) dada*, (1) 

Pi,j oc / (/i(a, a*) - /o,i) • {Ij{a,a*) - /qj) 

•/ a, a* 

X W{a,a*) dada*, (2) 
where a, a* are vacuum amplitudes of the (relevant set 



of) frequency modes [ig at the entrance of the crystal, 
W{a, a*) is obtained as the Wigner transform of the vac- 
uum state, Ii{a,a*) is the field intensity (for that mode) 
and /o,i the mean intensity due to the vacuum ampli- 
tudes, both at the entrance of the i-th detector (see .1^): 

/o = / Io{a,a*) W{a,a*) dada*. (3) 

J a.a* 

The last three expressions would in principle allow us to 
identify the vacuum amplitudes with a vector of hidden 
variables A € A (A is the space of events or probabilistic 
space), with an associated density function p(A), 

X = a,a*, piX) = W{a,a*). (4) 

Now, for instance in Q it is already acknowledged that 
([I])-© cannot, because of the possible negativity of the 
difference Ii{a, a*) — /o,i, be written as 

P^ ^ [ Pi{det\\) p{X) dX, (5) 

J A 

= / P^.,idet\X) p(A) dX, (6) 

J A 

where naturally Pi{det\X), Pij{det\X) should stay posi- 
tive (or zero) always. This last is our point of departure: 
in this paper we propose a reinterpretation of the for- 
mer marginal and joint detection probabilities based on 
a certain manipulation of expressions ([ij-lll]). 

The paper is organized as follows. In SeclUwe will con- 
sider a setup with just a source and two detectors; once 
that is understood, the interposition of other devices be- 
tween the source and the detectors poses no additional 
conceptual difficulty though it is nevertheless convenient 
to address it in some detail: this will be done in Sec. 
mil The calculations in these two sections find support 
on the proofs provided in Appendix [ll and stand for our 
main result in this paper. Sec. IIVI explores the ques- 
tion of "a-factorability" in the model, not only from the 
mathematical point of view but also providing some more 
physical insights on its implications. 

Up to that point the novel points of the paper are 
made and its results are self-contained; it is nevertheless 
natural to extend our analysis to some of their further 
implications (amongst these, the consequences for Bell 
tests of local-realism), in Sec. |Vl where we also include a 
preliminary approach to questions regarding the physics 
of the real detectors. Finally, overall conclusions are pre- 
sented in Sec. IVI[ and some supplementary material is 
provided in Appendix |21 which may not only help make 
the paper self-contained but perhaps also contribute to 
clarify some of the questions addressed, in particular the 
non-factorability issue. 

II. REINTERPRETING DETECTION 
PROBABILITIES 

Our reinterpretation of expressions involves 
two steps: 
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FIG. 1: Wigner-PDC scheme: photon pair generation, po- 
larizing beam sphtters (PBS) and detectors. Only relevant 
inputs of Zero-Point vacuum field (ZPF) are represented in 
the picture: "relevant" can be understood, in a classical wave- 
like approach, as "necessary to satisfy energy-momentum con- 
servation" for the (set of) frequency modes of interest; in a 
purely quantum electrodynamical one we would be talking 
about conservation of the commutation relations at the empty 
exit channels of the devices. Besides, those new ZPF compo- 
nents introduced at the empty exit channels of PBS's 1 and 
2 can alter the detection probability for the signal arriving to 
the detectors (for instance giving rise to its "enhancement"); 
that signal is on the other hand determined (amongst other 
hidden variables) by the ZPF components entering the crys- 
tal. Figure: courtesy of A. Casado. 



(A) Knowing that the following equality holds (from here 
on we drop detector indexes when unnecessary), for some 
real constant ("marginal"). 



K(m) / (^(«, a*) — /o) W{a,a*) dado* 

= [ P{det\X) p{X) dX, (7) 

J A 

we realize we do not need to assume 

P{det\X) = i^(„,)-(/(a,a*)-/o), (8) 

as a necessary, compulsory choice; it would be enough to 
find some f{x) > 0, satisfying 



^(m) / a*)-Io) W{a, a*) dada* 

f{I{a,a*)) W{a,a*) dada* 



so we can then safely identify 

P{det\X) = f{I{a,a*)), 



(9) 



(10) 



with f{I{a,a*)) > 0, V/(a,a*). This last is nothing 
but solving a linear system with only one restriction and 
an infinite number of free parameters, whose subspace 



of solutions intersects the region /(/(a, a*)) > Va, a*: 
see Appendix [TJ 

Of course, once at this point, joint detection probabil- 
ities would in principle impose additional restrictions on 
f{x); for detectors i,j we would have, for some constant 
("joint"). 



■^0) / (-^*("' "*) - Io,i) ■ a*) - Io,j) 

X W{a,a*) dada* 
f{Ii{a,a*)) ■ f{Ij{a,a*)) W{a,a*) dada*. 

(11) 



(B) In principle, it is unclear whether all possible /(a;)'s 
satisfying ([9|) can be made to meet pT|) . therefore we 
propose: instead of demanding from the beginning, 
we define a new function r(a;, y), so that 



^U) I - Io,i) ■ {Ij{a,a*) - Io,j) 

X W{a,a*) dada* 
r{Ii{a,a*),Ij{a,a*)) W{a,a*) dada*, 



(12) 

so we can now identify 

P^,Jidet\X) = r(/,(a,a*),/j(a,a*)). (13) 

In Appendix Q] we have shown that it is always possible 
to find some suitable < r(a;, y) < 1, for all x = Ii{a, a*) 
and y = Ij{a, a*), i.e., for all possible pairs a, a* . 

The absence of factorability on the a's, i.e., the fact 
tliat (we may have) 

T{I,{a,a*),Ij{a,a*)) ^ 

/(/,(«,«*))■/(/, (a, a*)), (14) 

may come as a surprise to some, given than Clauser- 
Horne factorability [19] on the hidden variable A is usu- 
ally taken for granted; this is a mistake We have 
tried to clarify the whole issue in Appendix!^ in general, 
we have shown that that non-factorability, aside from 
perfectly legitimate from the mathematical perspective, 
is in this case possibly related with the need to include 
new hidden variables in the model. ■ 

It is now convenient, for the following sections, to make 
our notation lighter, with a = a, a*, and also to define 

/,(«) = /(/.(a)), (15) 
f,,,(a) ^ mia),I,{a)). (16) 



III. ADDING INTERMEDIATE DEVICES: 
POLARIZERS, PBS'S... 

Once we place one or more devices between the crystal 
and the detectors, typically polarizers, polarizing beam 
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splitters (PBS) or other devices to allow polarization 
measurements (such as in Fig. [IJ, in general we cannot 
any longer describe the fields between both with only one 
set {a} of mode-amplitudes; we need to redefine our a's 
as now associated to a particular position r (they do not 
any longer determine a frequency mode for all space [l8l|). 
Hence, we will now have 



a(r) = {ak,^(r),ak,^(r)}. 



(17) 



and, letting be the position of the source (the crystal), 
and r; the position of the i-th polarizer or PBS (or any 
other intermediate device), we will also redefine 



(18) 



with ai including the relevant amplitudes at the (empty) 
exit channels of that i-th. intermediate device. With 
as,ai corresponding each to (a set of) modes with dif- 
ferent (sets of relevant) wavevectors {ks} and {k^} plj . 
we can then regard them as two (sets of) statistically 
independent random variables, with all generality. The 
intensity at the entrance of detector i-th will therefore 
depend now not only on but also on a^: 

Pi oc {Ii{as,ai) - Io,i) 

{h{as,at) - Iqa) W{as) W{at) dasdai. 

(19) 

Once more, something like that can always be rewrit- 
ten (see former section) , for some suitable and positively 
defined f'{x), as 



P, 



f'{Ii{as,ai)) W{as) W{a.i) dasdai 



Integrating on ai we would obtain 

P, = / f[{as) W{as) das 



(20) 



(21) 



from where we define a new function fl{as). On the other 
hand, for joint detections we would have 



Pi, J oc 



{Ii{as,ai) ~ Iqa) ■ {Ij{as,aj) - Iqj) 

X W{as) W{ai) W{aj) dagdaidaj, 

(22) 



which again can always be rewritten (again see former 
section), for some positively defined T'(x,y), as 



P, 



r'{I,{as,ai),Ij{as,aj)) 

X W{as) W{at) W{aj) dagdaidaj, 

(23) 



and integrating on ai, aj we would obtain 

P,,j = I Tl,{as) W{as) das, (24) 

J as 

from where we can again define yet another new proba- 
bility density function T^j{as). 

It is interesting for the sake of clarity to compare the 
two situations: with (primed functions) and without po- 
larizers (unprimed). It is easy to see that, because the 
detector only sees the intensity at its entrance channel, 
clearly (we drop the "s" subscript for simplicity). 



fix) = fix), Wx, 
r'(a;,y) = T{x,y), \/x,y. 



(25) 
(26) 



while, naturally, in general //(a) ^ fi{oi), as well as 
r^j(a) ^ rij(a)- 

IV. ON NON-FACTORABILITY 

A. Mathematical analysis 

Let us go back to the case with just the source and the 
detectors; we will soon see the following does nevertheless 
also apply when polarizers or other devices are added to 
the setup, just the same. According to our reasonings in 
App. [2l and using ([15]) -([16]), we now realize that there 
is no way to avoid 



r»j(a) = h{a) ■ fj{a), 



(27) 



unless we introduce some additional dependence of the 
kind /(a) f{a,^), so that then 



(28) 



where we add a second "hat" to avoid an abuse of nota- 
tion, and where fj, stands for a new set of random vari- 
ables. This f{a,^) should be interpreted as a detection 
probability conditioned to the new vector of random vari- 
ables /J,, i.e. 



f{a,n) = P{det\a,^i). 



(29) 



We will impose further demands on f{a,fi), defining 



f{a,fi) = f{I{a,n)), 



(30) 



something forced by strictly physical arguments: the 
choice f(I{a,fj,)) must prevail over other possible ones 
- for instance f{I{a), ^) - due to the need to respect the 
dependence of the probabilities of detection (conditioned 
to a or not) alone on the intensity that arrives to the de- 
tector, and nothing else. Now, with the density function 
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we could write 

P{det\a) ~ / P{det\a, ij) p^{iJ,) dfi, (31) 

Pij{det\a) = / Pi,j{det\a, fi) Pfj,{p) dp,, (32) 
J 

allowing us to recover our former definitions (jl5|) ~ (|16p : 

f{a) = P{det\a), (33) 
f,,,(a) =P,,,(det|a). (34) 

For joint detections, the additional variable p is particu- 
larly relevant because, we will always have that while 

Pij{det\a, p) = Pi{det\a, p) ■ Pj{det\a, p), (35) 

in general 

P^j{det\a) ^ P,{det\a) ■ Pj(det\a), (36) 
or we could equivalently say that while necessarily 

t^J{a,p) = fi{a,p) ■ fj{a,p), (37) 

in general 

f, (a) ^ /,(«)• /,(«), (38) 

where of course 

f*j(a) = / Ma,p) ■ fj{a,p) p^{p) dp. (39) 

To conclude this section, we recover the case with in- 
termediate devices; due to a^, aj being, as defined, inde- 
pendent from one another and also from as, our hypo- 
thetical "flag" p cannot be associated with none of them. 
Therefore, in general, and in principle, not only 

f^,^.(a)^ /;(«)• /;(«), (40) 

but also Ti^j{a) ^ fi{ot) ■ fji<^) either. We use "in prin- 
ciple" because this question is not yet analyzed in de- 
tail; we now see crystal clear, though, that this possible 
non-factorability on a's is nothing more than an internal 
feature of the model's mathematical structure, bearing 
no relevance in regard to its double-sided compatibility 
(or absence of it) both with local-realism and the quan- 
tum predicted correlations. Indeed, the crucial element 
in that sense is the intervention of new random vacuum 
amplitudes introduced by the new devices in the setup, 
as we will see in Sec. El 

B. A physical interpretation 



Pi{det\a) ■ Pj{det\a), where p would be a new, additional 
HV (or vector of HV's) and should, in the ultimate term, 
correspond to some physically relevant variable that we 
have so far overlooked, perhaps containing additional in- 
formation about the state of the crystal and the laser. 
Indeed, if it was eventually necessary to introduce a new 
variable p, physical intuition suggests that it could prob- 
ably be associated to the phase of the laser complex am- 
plitude [2^; with the laser in a coherent state, it would be 
justified to approximate p^i{p) ~ S{po), where po would 
be the phase or the quasi-classical wave. Perhaps it is 
needless to add that in such a case expressions such as 
(p9| or a similar one for r"s would still be fully opera- 
tive, as long as the dependence on po is "absorbed" in 
the functional dependence at the left side of the equation. 
Anyway, whether this p has an obvious physical mean- 
ing or not should not bother us too much here, as it does 
not affect at all the correctness and/or the implications 
of our mathematical results. 



V. COMPLEMENTARY QUESTIONS 

A. Wigner-PDC's local realism vs. quantum 
correlations 

Though former mathematical developments are fully 
meaningful and self-contained on their own, yet it would 
be convenient to give some hints on how a local-realist 
(LR) model can account for typically quantum correla- 
tions, which are known to defy that very same local re- 
alism (LR). In the first place and as a general answer, 
what the Wigner-PDC picture proves is that LR is re- 
spected by a certain subset of all the possible quantum 
states, specifically the ones that can be generated from a 
non-linear mix of the QED- vacuum (which therefore acts 
as an "input" for the model) with a quasi-classical (a 
high-intensity coherent state), highly directional signal, 
the laser "pump" (which indeed enters in the model as 
a non-quantized, external potential @). Moreover, such 
a restriction is clearly not arbitrary at all, since it arises 
from a very simple quantum electrodynamical model of 
the process of generation of polarization-entangled pairs 
of photons from Parametric Down Conversion (PDC): see 
for instance eq. (4.2) in 



B. Detection rates and "efficiencies" 

Aside from subscripts, we will now also drop "hats" 
and "primes" for simplicity; of course the fact that 
Tij{a) may not be in general a-factorisable, 



Recapitulating, we have that while Pij{det\a, p) = 
Pi{det\a, p) ■ Pj{det\a,p), however Pij{det\a) ^ 



r»j(a) / Ma) ■ fj{a), 



(41) 
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does not at all mean that it cannot well satisfy 



r,.,(a) W{a) 



da 



fi{a) W{a) da 



f.j{a') W{a') da' 



(42) 



i.e. (let us from now use superscripts "W" and "exp" to 
denote, respectively, theoretical and experimental detec- 
tion rates): 



P^J\det)^P^'^\det).P^'^\detl 



(43) 



which is indeed the sense in which the hypothesis of "er- 
ror independence" is introduced, to our knowledge, in 
every work on LHV models [25| . This sort of conditions 
over "average" probabilities (average in the sense that 
they are integrated in the hidden variable, may that be 
a alone or also some other one) are the only ones that 
can be tested in the actual experiment; there, we can 
just rely on the number of counts registered on a certain 
time- window AT, and the corresponding estimates of the 
type 



Pt^P\det) 



n. 



joint det. in AT 
, marg. det. (j) in AT 



(44) 



Now, if we wish to include some additional uncertainty 
element reflecting the technological limitations (a "de- 
tection efficiency" parameter), what we have to do is to 
redefine the overall detection probabilities as 



i 

r,{exp) 



(det) 
(det) 



(det), 



as well as 



i 

,{exp) 



{det\a) = f], ■ Pl'^\det\a) 



{det\a) 



(45) 
(46) 

(47) 
(48) 



where < ?}i < 1 would play the role of such an (al- 
leged) detection efficiency, the "hats" remarking the fact 
that the customary definition of the analogous quan- 
tity in QInf involves not only our 7)'s but also the non- 
technological contribution. From the point of view of 
the experimenter it is very difficult to isolate both com- 
ponents (Glauber's theory [l5j does not predict a unit 
detection probability even for high intensity signals); we 
should perhaps then confine ourselves to the term "ob- 
servable detection rate" instead of using the clearly mis- 
leading one of "detector inefficiency" . 



C. Consequences on Bell tests, their 
supplementary assumptions and critical efficiencies 

That said and going to a lowest level of detail, states 
in such an (LR) subset of QED can still indeed exhibit 



correlations of the class that is believed to collide with 
LR, yet the procedure through which they are extracted 
from the experimental set of data does not meet one of 
the basic assumptions required by every test of a Bell 
inequality: they do not keep statistical significance with 
respect to the physical set of "states" or hidden instruc- 
tions '25']). To guarantee that statistical significance we 
must introduce some of the following two hypothesis: 

(i) all coincidence detection probabilities are independent 
of the polarizers' orientations 0i, 4)j (this is what we call 
"fair-sampling" [26|, for a test of an homogeneous in- 
equality i27j), which implies 



Pij{det\(j)i,4}j,a) = Pij{det\as), 



(49) 



where of course (see Sec. IIII|) a = as (B ai (B aj, and 
where we recall that , would determine which 
vacuum modes amongst the sets , aj would intervene 
in the detection process, 



(ii) the interposition of an element between the source 
and the detector cannot in any case enhance the probabil- 
ity of detection (the "no-enhancement" hy pot hesis [l^l, 
needed to test the Clauser-Horne ineaualitv [1911 , and pre- 
sumably every other inhomogeneous one [29|), 



Pi{det\(f>i,a) < Pi{det\oo,as 



(50) 



with 00 denoting the absence of polarizer and with a = 
as ® ai this time. 

Following our developments in Sec. IIIII one can easily 
see that (i) is not in general true, and according to |13j 
neither is (ii). I.e., whenever states of light are prepared 
so as to produce the sort of quantum correlations that are 
known to defy LR, these last come supplemented with the 
necessary features that prevent it from happening... how 
could it not? 

Yet, a mere breach of (i) or (ii) is not enough to assert 
the existence of a Local Hidden Variables (LHV) model, 
which is an equivalent way of saying that the results of 
the experiment respect LR: it is more than well known 
that this can only happen for certain values of the ob- 
served detection rates jS^H^]. However, from the point 
of view of (|T7|) - (P5)) . and given the fact that, as proven in 
Sees nil and IHIl the Wigner-PDC is in all circumstances 
in accordance with expressions (|54)) - (|55|) . and hence to 
all possible Bell inequalities whether our 77 's are equal or 
less than unity, the so-called "critical efficiencies" would 
merely stand, at least as far as PDC-generated photons 
are concerned, for bounds on the detection rates that 
we can experimentally observe (these last in turn con- 
strained by the only subset of quantum states that we 
can physically prepare). 
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D. An approach to realistic detectors and average 
ZPF subtraction 



In this section we make a slight approach on some phys- 
ical considerations with the aim of showing that there is 
plenty of room for a suitable physical interpretation of 
the model, even when that is not strictly necessary for 
the coherence of our results, at least from the purely 
mathematical point of view. Indeed, we have already 
descended to the physical level when we established, in 
former sections, the dependence of our /'s and F's solely 
on the intensity arriving at the detector. 

Expectable behavior for a physical device would typ- 
ically include a "dead-zone", an approximately linear 
range and a "saturation" at high intensities (this is in- 
deed the kind of behavior suggested for instance in [IB]); 
amongst other restrictions this would imply, for instance, 
/(/(a)) = 0, when I{a) < Iq (this last a threshold that 
may even surpass the expectation value of the ZPF in- 
tensity), as well as T{Ii{a),Ij{a)) = either for Ii{a) 
or Ij{a) below Iq. Neither these restrictions nor other 
similar ones would in principle invalidate our proofs in 
Appendix [1] which seem to provide room enough to sim- 
ulate a wide range of possible behaviors; however, we 
must remark that none of our /'s and F's can ultimately 
be considered as fully physical models, due to the fact 
that they represent point-like detectors (the implications 
of such an over-simplification may become clearer in the 
light of our following point here). 

In close relation to the former, we also propose here a 
simple physical interpretation of the term —Iq appearing 
in the expressions for the detection probabilities. From 
the mathematical point of view such subtraction arises 
from a mere manipulation of Glauber's original expres- 
sion [15l] . From the physical one, a realistic interpretation 
would be more than desirable, as that subtraction of ZPF 
intensity is for instance crucial to explain the absence of 
an observable contribution of the vacuum field on the de- 
tectors' rates [3l|; of course we mean "explain in physical 
terms" ; from the mathematical point of view our model 
here already predicts a vanishing detection probability 
for the ZPF alone. Our suggestion is that the —Iq term 
must be (at least) related to the average flux of energy 
going through the surface of the detector in the opposite 
direction to the signal (therefore leaving the detector). 
That interpretation fits the picture of a detector as a 
physical system producing a signal that depends (with 
more or less proportionality on some range) on the total 
energy (intensity times surface times time) that it accu- 
mulates. For the moment we just say "at least related" , 
bearing in mind that to establish such association we 
would first have to refine the point-like model of a de- 
tector which stems directly from the original Glauber's 
expression [isj . 



VI. CONCLUSIONS 

We have shown that the Wigner description of 
PDC-generated (Parametric Down Conversion) photon- 
entanglement, so enthusiastically developed in the late 
nineties but then ignored in recent years, can be for- 
mulated as entirely local-realistic (LR) . A formalism that 
is one-to-one with a quantum (field-theoretical) model of 
the experimental setup can be cast, thanks to an addi- 
tional manipulation (also one-to-one), into a form that 
respects all axiomatic laws of probability 32] , and there- 
fore LR, as defined for instance in Sec. 1 2 al The original 
quantum electrodynamical model takes as an input the 
vacuum state, which accepts a well defined probabilistic 
description through the so-called Wigner transform: this 
is the fact that the analogy with a local-realistic theory 
is conditioned to. What we call Wigner-PDC accounts, 
then, for a certain subset within the space of all possi- 
ble QED-states, determined by a particular set of ini- 
tial conditions and a certain Hamiltonian governing the 
time-evolution ^^\, restrictions that seem to guarantee 
(according to us here) the compatibility with LR. 

Aiming for the maximum generality, as well as to 
avoid some possible (still under examination) difficulties 
with factorisable expressions, we have renounced to what 
we call a-factorability of the joint detection expressions. 
Such a choice is not only perfectly legitimate [20^, but 
may also be supported by a well feasible interpretation 
(see Appendix IIV[) : nevertheless, further implications of 
that non-factorability on a will also be left to be exam- 
ined elsewhere. Again, whatever they finally turn out to 
be, they are also irrelevant for our main result in this 
paper: the Wigner-PDC formalism can be cast into a 
form that respects all axiomatic laws of probability for 
space-like separated events. 

Neither does the explicit distinction between the cases 
where polarizers (or other devices placed between the 
source and the detectors) are or not included in the setup 
introduce any conceptual difference from the point of 
view of our main result; however, the question opens 
room to remark some of the main differences of the 
Wigner-PDC model (actually, also its Hilbert-space ana- 
logue) with the customary description used in the field 
of Quantum Information (QInf): here, each new device 
introduces noise, new vacuum amplitudes that fill each 
of its empty polarization channels at each of its exists, in 
contrast to the usual QInf "black box" , able to extract 
polarization information from a photon without any in- 
deterministic component. 

As a matter of fact, those additional random compo- 
nents hold the key to explain the variability of the detec- 
tion probability (see note 13]) that is necessary, from the 
point of view of Bell inequalities, to reconcile quantum 
predictions and LR (see Sec. IV Ap . In particular, the im- 
mediacy with which the phenomenon of "detection prob- 
ability enhancement" arises in the Wigner-PDC frame- 
work would suggest that this may be after all the right 
track to understand why after several decades the min- 
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imum detection rates (or, in QInf terminology, critical 
detection efficiencies) that would lead to obtain conclu- 
sive evidence of non-locality are yet to be reached [soj . 

In Sec. IVDI we have done a first, general approach on 
the question of whether our reinterpretation of the de- 
tection probabilities is consistent with the actual physi- 
cal behavior of detectors. A closer look to this question 
is nevertheless left out of the scope of the paper; former 
proposals [l^ in regard to this issue aimed perhaps too 
straightforwardly to the physical level, while they did 
not even guarantee consistency with the framework we 
have settled here. In close relation with the former, we 
have also suggested a possible simple physical interpre- 
tation of /o-subtraction taking place in the expressions 
for the detection probabilities: work in any of this two 
directions would anyway require a departure from the 
point-like model of a detector on which we have focused 
here. Summarizing, we have shown that a whole family 
of detection models 

Mdet = {/(/(a,^)),r(/,(a,/i),/j(a,M)}, (51) 

can be found, consistent both with the quantum mechan- 
ical expressions from the Wigner-PDC model and LR. 
A close examination of the constraints coming from the 
physical behavior of the detectors and other experimen- 
tally testable features is left as a necessary step for the 
future, with the aim of establishing a subset physically 
feasible ones; nevertheless, we have the guarantee that 
all of them produce suitable predictions, as so does their 
quantum electrodynamical counterpart. 

Yet, even at the purely theoretical level some other 
features remain open too: as a fundamental one, to 
what extent the model requires what we have called non- 
factorability on a's. Anyway and as a final conclusion, 
though our work here does not imply that we have al- 
ready achieved an entirely proficient description of PDC- 
based photon entanglement, at least for us is already ev- 
ident that, up to the point we have come up to, that 
description is indeed completely local-realistic. 

It is more than proven, yes, that QM is sometimes in- 
compatible with local- realism, but... is really every quan- 
tum state physically realizable? Does QM describe na- 
ture or is it just a subset of it (the one that is compatible 
with local-realism) that really does? After all, QM is just 
a theory, a theory that provides a formalism upon which 
to build models, models than can (and should) be refined 
based on experimental evidence; quantum entanglement 
seems to manifest in many of its "reasonable" features 
but, at least as PDC is concerned and so far, whenever 
local-realism would seem to be challenged new phenom- 
ena appear preventing it. Such are "unfair sampling", 
"enhancement" (as a particular case of variable detec- 
tion probability) and over all detection rates low enough 
to open room for the former two. To explore, and ex- 
haust if that is the case, alternative routes such as the 
one here is not only sensible but also necessary. 
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Appendix 
1. Auxiliary proofs 

Lemma 1: There always exists fix) satisfying with 
< f{x) < 1, for X = /(a, a*) and all possible pairs 
a, a* . 

Proof: We have a linear system with one restriction and 
infinite independent variables: N"^ (N — > oo) real num- 
bers f{I{a, a*)), one per each pair a, a* . Coefficients are 
given by the W{a, a*)'s, and the independent term is the 
left-side term of -Ps > 0. We consider the real, vecto- 
rial space V where a particular solution of the problem is 
represented associating one /(/(a, a*)) per coordinate: 
dim(V) = (again with N — ^ oo). If they do exist, 
compatible solutions /(/(a, a*)) for the (unique) restric- 
tion will conform a linear manifold A4 C V. Now we have 
to see: 

(i) Due to W{a,a*) > for all pairs a, a* Q, and 
Ps > too, A4 cannot be parallel with any of the coor- 
dinate hyper-planes in V; i.e., Ai intersects all of them, 
defined each (each one for a pair a, a*) by the equation 
/(/(a,a*)) = 0. 

(ii) Moreover, for > 0, has a non-trivial intersec- 
tion (more than one point, the origin) with all coordinate 
planes inside the first hyper-quadrant Viq (the subregion 
Viq C V given by restricting V to f{I{a,a*)) > 0). This 
can be seen, for instance, determining the point of cross- 
ing with the axes: doing all f's zero except one (always 
possible because all the Ws are strictly above zero), we 
can see the crossings always take place at the positive 
half of the corresponding coordinate axis, therefore at 
the boundary of Viq. For Pg = 0, the solution for the 
system is trivial. 

With (i) and (ii), it is clear that under the restriction 
(under the set of inequalities) f{I{a,a*)) < 1 Va, a*, 
the set of admissible solutions is still not empty. This is 
guaranteed with 




W{a,a*) dada* = 1. (52) 



and < 1- Here we give an inductive reasoning; think 
in 3 dimensions: what we have is ax + by + cz — d, with 
a,b,c > and d > 0; the coordinates of the point of 
intersection with x — y — z are given (obviously all) by 
the quantity c?/(a -f 6 + c). The extension to infinite di- 
mensions, topological abnormalities all absent, is direct. 
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with a,b,c = W's, d = Ps and a + b + c = 1, and finally 
-Ps < 1, guaranteeing the existence of valid solutions. ■ 

Lemma 2: There always exists T{x,y) satisfying Iil2\} . 
with < T{x,y) < 1, for x — Ii{a,a*), y — Ij{a,a*) 
and all possible pairs a, a* . 

Proof: formally identical with Lemma 1. ■ 
Lemmas 1 and 2 are directly applicable, aside from Sec. 

in to Sec. EH 



2. Basic concepts revisited 

We first briefly revisit the concepts of locality, deter- 
minism and factor ability; a good understanding on these 
concepts is crucial for the main results of the paper, what 
makes this review not only convenient but almost un- 
avoidable, especially given the presence of some confusion 
in the literature. In any case, it shall be clearly under- 
stood that Clauser and Home's factorability flS] is not 
a requisite for local-realism. 



a. Locality and realism 

A theory predicting the results of two measurements 
A and B that take place under causal disconnection (rel- 
ativistic space-like separation) can be defined as local if 
and only if we can write 



variable (HV) A (summarizing the state of that system) , 
if (and only if) 

P(M = TOfc|A) e {0,1}, Vfc,A, (56) 

which allows us to write 

M = Af(A), (57) 

and indeterministic iff, for some A, some k' , 

P{M = mk'\X)^ {0,1}, (58) 

i.e., at least for some (at least two) of the results for at 
least one (physically meaningful) value of A. 

Now, indeterminism can be turn into determinism, i.e., 
(1551) can into ((5E)) . by defining a new hidden variable /x, 
so that now, with 7 = A © /x: 

P(Af = mfe|7) e {0,1}, Vfc,7, (59) 

which means we can write, 

M = M(A,Ai), (60) 

a proof that such a new hidden variable fi can always be 
found (or built) given in [35|. 

So far, then, our determinism and indeterminism are 
conceptually equivalent, though of course they may cor- 
respond to different physical situations, depending for in- 
stance on whether 7 is experimentally accessible or not. 



A^A{X,c^a), B = B{\,<Pb), (53) 

where A is a (set of) hidden (or explicit) variables defined 
inside the intersection of both light cones, and (jjAjf^B 
are another two other sets of variables (amongst them 
the configurable parameters of the measuring devices) 
defined locally at A and B, respectively, and causally 
disconnected from each other, i.e., 

P{A^a\X,(bA,^B)^ P{A^a\X,(bA), (54) 
PiB - b\X, (t)A, 0b) = P{B = 6|A, 0b). (55) 

These last two expressions are usually taken as a defini- 
tion of local causality (33 |. 

Now, realism simply stands for A (and (pA, i'B as well) 
having a well defined probability distribution. A set of 
physical observables corresponding to a particular quan- 
tum state can be sometimes described by a well defined 
joint probability density (such is the case of field ampli- 
tudes in any point of space for the vacuum state in QED); 
for other quantum states that is not possible though. 



c. Factorability 

Let now M = {Mi} be a set of possible measurements, 
each with a set {rn^./t} of possible outcomes, not neces- 
sarily isolated from each other by a space-like interval. 
We will introduce the following distinction: we will say 
M is 

(a) X-factorisable, iff we can find a set {^i} of random 
variables, independent from each other and from X too, 
such that 

A* = 0C., (61) 

i 

and ([55]) holds again for each Mi on 7^ = A © : 

P(M, =m,,fc|7,) G {0,1}, Vz, Vfc,7„ (62) 

this last expression meaning of course that we can write, 
again for any of the M^'s, 

= M,(A,C»). (63) 



b. Determinism 

A measurement M upon a certain physical system, 
with k possible outcomes ruk, is deterministic on a hidden 



(b) non X-factorisable, iff (|62p is not possible for any set 
of statistically independent ^^'s. 

We will restrict, for simplicity, our reasonings to just 
two possible measurements A, B A4, with two possible 
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outcomes, A, B € {+1,-1}, all without loss of gener- 
ality. We have seen that, as the more general formula- 
tion, we can always write something like A — A{X,£^a), 
B = B{\^b). 

Lemma 1~ 

(i) If A and B are deterministic on A, i.e., (j56p holds for 
A and B, then they are also A-factorisable, i.e., 

P{A ^a,B^ b\X) = P{A = a\X) ■ P{B = 6|A), (64) 

for any a,b ^ {-|-1, — 1}. Eq. ([M)) is nothing but the 
so-called Clauser-Horne factorability condition [1SJ|. 

(ii) If A and B are indeterministic on A, i.e., if (|56|) does 
not hold for A, then: for some /x (always possible to find 
[ssj ]! such that now (f59| holds for 7 = A © ^, A, i? are 
7-factorisable, 

P{A = a,B ^ &I7) = P{A ^ a|7) • P{B = 6I7), (65) 

i.e., ((64|) holds for 7, this time not necessarily for A. 

(iii) Let ([62|) hold for A,B, on X,^a,^b- if ^,£,a,£,b are 
statistically independent, (hence, A and B are what we 
have called A-factorisable), then (|64)l holds for A, not 
necessarily on the contrary. 

Proof- 

(i) When ([56]) holds, for any A and any a,b E {+1, —1}, 
P{A = a\X),P{B = &|A) e {0,1}, from where we can, 
trivially, get to (|64l)- ■ 



(ii) It is also trivial that, if ([5^ holds, (IMl) can be recov- 
ered for 7. That the same is not necessary for A can be 
seen with the following counterexample: suppose, for in- 
stance, that for A = Ao, either A = B = 1 ot A — B = —1 
with equal probability. It is easy to see that 

P(A^B^ l\Xo) ^ P{A = l|Ao) • P{B = l|Ao), (66) 



numerically: \ ^ \- ■ 

(iii) We have, from independence of A,^yi,^B, and work- 
ing with probability densities p's: P\{X,£_Ai£,b) — P\W ■ 
Pa{^a) • P 3(1,3), which we can use to write 



P{A = a,B = b\X) = J P{A = a,B = b\X, Cs) 

X pa{S.a) ■ Pb{£,b) d^Ad£,B- 
(67) 



and now with the fact that we can recover ([5^ for A (B) 

owfA^ XS)S,A (75 = A © Cb), 



P{A = a,B = b\X) 

= J P{A^a\X,U)-P{B = b\X,^B) 

X Pa{£,a) ■ Pb{£,b) d£,Ad^B 
= J P{A = a\X,U)-PA{U) dU 

X J P{B^b\X,^B)-PB{^B) d^B 



P{A ^ a\X) ■ P{B = 6|A). 



(68) 



On the other hand, let X^S^a,£,b be not statistically in- 
dependent: we can set for instance, as a particular case, 
= /i, Vi, therefore reducing our case to that of (|59l) . 
Once this is done, our previous counterexample in (ii) is 
also valid to show that factorability is not necessary for 
A here. ■ 
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value of the number of photons is zero and therefore such 
rate vanishes: experimentally observed dark counts are 
usually associated to thermal effects or electronic noise. 
Besides, a saturation effect is introduced, which on the 
other hand is a feature of all real detectors, but in this 
case it appears at much lower intensities (for energies 
lower than one photon). Consistency with experimental 
data demands compliance with certain relations provided 
in that same paper [l6l ]. 

Indeed, the laser is described by a coherent state, which 
for a high intensity signal means that it can be regarded 
as quasi-classical wave with well defined amplitude and 
phase (and introduced as a non-quantized external po- 
tential, see 0] for instance). 

The phase of this complex amplitude is for instance rele- 
vant in expressions (4.10) from 3|, and has the potential 
to interfere constructive or destructively with the a's, in- 
creasing or decreasing the overall intensity of the signal, 
therefore modifying the detection probabilities. From this 
point of view, such phase cannot be at all excluded from 
the vector of relevant hidden variables and therefore non- 
factorability (of detection probabilities) in a is the most 
natural feature to expect (though not strictly necessary), 
see our detailed analysis in Sec. IIVI 
R. Risco-Delgado, private communication. 
For reference on how, and under which hypothesis LHV 
models are built: 

(i) A. Cabello, D. Rodriguez, I. Villanueva. Phys. Rev. 
Lett. 101, 120402 (2008), 

(ii) A. Cabello, J. -A. Larsson, D. Rodriguez. Phys. 
Rev. A 79, 062109 (2009). 

The fair sampling assumption is already formulated in 
[37I ]: "given a pair of photons emerges from the polariz- 
ers, the probability of their joint detection is independent 
of the polarizer orientations" . 

We believe the distinction between homogeneous and 
mhomogeneous inequalities was first introduced by M. 
Horne himself {ref), then later updated by E. Santos 



(ref). We will supersede previous definitions with our 
own here: we will call mhomogeneous inequalities all the 
ones involving coincidence rates of different order (for 
instance, marginal and joint ones), homogeneous on the 
contrary. 

Of course, if for instance a certain subset of the observ- 
ables involved in the inequality does not have any un- 
certainty associated to its detection, the inclusion of in- 
homogeneous terms may not lead to the requirement of 
supplementary assumptions; this is not the case of exper- 
iments with photons, anyway. 

[28] The no-enhancement assumption [l^: "for every emis- 
sion A, the probability of a count with a polarizer in place 
is less or equal to the probability with the polarizer re- 
moved", where A is the (hypothetical) hidden variable 
expressing the state (at least the initial one) of the pair 
of particles. We remark: "for every emission A". 

[29] An inhomogeneous Bell inequality requires the es- 
timation of coincidence rates of different order (for in- 
stance single and double, or double and triple); due to 
the fact that there is no way to identify if the whole set 
of particles have been simultaneously emitted and then 
gone undetected, or they were simply never emitted at 
all, the test would always require supplementary assump- 
tions; of course from a wave-like perspective, the issue is 
even more evident. In any case, marginal detection rates 
would not be experimentally accessible unless we assume 
"independent errors" at the detector for every and each 
of the "states" in the model (i.e., detection probabilities 
cannot depend on the hidden variables: a's, other such 
as the phase of the pump. etc). 

[30] Let us consider a Bell experiment and let r; what people 
in QInf define as "detection efficiency" [25| (from our 
point of view this is clearly a misleading term, we should 
call it simply "detection rate" , whether is reduced or not 
by technological imperfections, see Sec. IV A")) . A possible 
Local Hidden Variables (LHV) model for the experiment 
is then composed by a set of states {s} with probabilistic 
weights (or frequencies) ps'- all restrictions the model 
must satisfy (either regarding the behavior of detectors 
or the quantum mechanical predictions) are linear in {ps} 
and the trivial solution rj = (all ps = but for the 
one that predetermines no detections at all) is always 
admissible, implying that there is always some T^crit > 
such that for all r] < r/crit the desired LHV can be built. 
The reasoning applies whether if the LHV simulating just 
simulates a violation of a particular Bell inequality or also 
every other quantum prediction for a given state and set 
of observables (though the addition of more restrictions 
may lower r/crit, as expectable). 

[31] Though from the mathematical point of view our model 
here already predicts a vanishing detection probability 
for the ZPF alone, a physical interpretation of the 
absence of an observable detection rate as a result of 
the vacuum fiuctuations, or at least that of a significant 
one, is still an open question. Following Santos's work, 
the absence of the ZPF from the observational spectrum 
could be justified, for realistic detectors, on the combi- 
nation of the following properties: 

(i) the already commented subtraction of the average 
ZPF intensity as a result of Glauber's expression in the 
normal order [l5| , 

(ii) spatial and temporal integration, involving the 
autocorrelation properties of the ZPF, 
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(iii) a "low band pass" frequency response. 
See, for instance: Emilio Santos, "How pho- 
ton detectors remove vacuum fluctuations", 
http://arxiv.org/abs/quant-ph/0207073v2 (2008). 
We remark again that this issue is in any case irrelevant 
for our main results here: it concerns models built 
outside of the quantum framework (which is not our 
case). Nevertheless, we conjecture that our proposed 
interpretation of (i) in Sec. IVII may be of use to make 
progress on the question. 
[32] Let A be a space of events, all possible probability assig- 
nations within the formalism, F{-) : A £ A — >■ F{X), 
satisfy < F{X) < 1 VA and, for any partition {A^} of 
A (assuming of course additivity on disjoint subsets of 
events), J2^F{A,) = 1. 

All Bell inequalities (at least those written in terms of 
probabilities) can be obtained as a derivation from these 
laws, plus some very simple supplementary assumptions: 
statistical independence of variables defined at distant 
locations. Of course. Bell inequalities written in terms of 
correlations can be also written in terms of probabilities. 
Many do no even need such supplementary assumptions 
(space- like separation ones or of another kind): such is 
the case of both the CH 19] and the CHSH ;37] inequal- 
ities. 

To illustrate this see for instance [3^'s derivation of the 
CH inequality; the CHSH inequality can on the other 
hand also be interpreted as a mere algebraic inequality 
on whatever four quantities taking values ±1. 



[33] W{a,a*) is a Gaussian 0|; hence W{a,a*) > Va,a*. 

[34] For a recent work reviewing these concepts see for in- 
stance: T. Norsen, "J.S. Bell's Concept of Local Causal- 
ity", http://arxiv.org/abs/0707.0401. 

[35] A possible (not unique) procedure to build /i is this: for 
each Mi, we define a new random variable ai and assign 
values for each pair {X,mi.k): (Ti = ai{X,mi^k), and now 
simply do = ®-(Ji. As built, tJi's are not necessarily 
independent from one another, nor are they necessarily 
independent from A. 

[36] A. Casado, private communication. 

[37] J.F. Clauser, M.A. Home, A. Shimony, R.A. Hoh. Phys. 

Rev. Lett. 23, 880 (1969). 
[38] J.-A. Larsson and J. Semitecolos, Phys. Rev. A 63, 

022117 (2001). 

[39] R. Risco-Delgado. PHD Thesis. Universidad de Sevilla 
(1997). 

[40] The PDC model is by construction restricted to a cer- 
tain subset of QED-states, obtained directly from a mix 
of the vacuum state (hence one with positive Wigner 
function) with a quasi-classical signal (the laser), and 
a time-evolution governed by a quadratic Hamiltonian 
(hence one that preserves the positivity of the Wigner 
function); for this last point, see for instance A. Casado, 
PHD Thesis: "Estudio de los experimentos de conversion 
parametrica a la baja con el formalismo de la funcion de 
Wigner" (Universidad de Sevilla, 1998). 



